Formules efficaces pour ¢(n) 


par 
Simon Plouffe 
27 novembre 2023 


Résumé 


Dans les Ramanujan Notebooks, on peut apercevoir une formule originale. 
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Qui m’a inspiré pour en trouver d’autres semblables bien aprés avec les puissances de z et ((n), n 
> 1. Je présente ici une nouvelle gamme de formules duméme genre mais nettement plus efficaces 
pour le calcul. Elles s’appliquent a toutes les valeurs de n. Ces nouvelles formules permettent 
d’évaluer les toutes les valeurs de la fonction Zeta avec une efficacité de 3.85 et 5.46 décimales par 
terme. En exploitant certaines de ces séries on obtient une formule explicite pour la n’iéme 
décimale de ¢(3). 


Summary 


In the Ramanujan Notebooks, we can see an original formula. 
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Which inspired me to find similar ones much later with powers of 1 and ¢(n),n > 1. Herel present 
a new range of formulas of the same kind, but considerably more efficient for calculation. They 
apply to all values of n. These new formulas allow all Zeta function values to be evaluated with an 
efficiency of 3.85 and 5.46 decimal places per term. 


Introduction 


Nous revenons sur ces formules trouvées en 2006, 
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On remarque un motif dans les exposants (le terme avec e™"), le motif est toujours (1, 2, 
4) mais ce n’est pas le seul. Dans le cas de ¢(5) le motif est plutét (4, 5, 10, 20). En fait en 
y regardant de prés, le motif (1, 2, 4) sont les diviseurs de 4 et (4, 5, 10, 20) est un sous- 
ensemble des diviseurs de 20. De plus, la série pour ¢(5) converge relativement 
rapidement puisque chaque terme donne 5.46 décimales. Le motif (4, 5, 10, 20) est par 


ailleurs répété avec tous les ((4n + 1). 


En effet, les coefficients pour n=5, 9, 13, 17, ... avec le motif (4, 5, 10, 20) persiste. 
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En reprenant l’une des formules de L. Vepstas [11], on peut trouver d’autres expressions, 
il en existe une infinité mais la plus efficace qui a été trouvée pour ¢ (3) est ; 
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Qui donne bien 3.86 décimales par terme. Elle est valide aussi aux valeurs ¢(4n — 1), les 
coefficients étant différents. 


Quant a la formule avec V7 et ¢(5), de la méme facon, s’applique aux valeurs ((4n + 1). 
Donc, la meilleure trouvée pour ¢(5) et les ¢(4n + 1) permet tout de méme d’avoir 5.46 
décimales par terme. L’efficacité augmetant avec l’exposant de n. A ma connaissance ces 
formules sont les plus simples et les plus efficaces connues méme si on doit évaluer la 
valeur de 7. 


La formule avec ¢(3) ci-haut est moins efficace que les formules de Amdeberhan et 
Zeilberger mais elle est plus simple a programmer. 


Une autre modéle générateur de formules est le suivant. 
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Ou A, B et C sont rationnels, par exemple avec ¢(3),¢(5) et ((7) ona 
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Elle est valide Vm > 3. 


Donc le 2&™¢ terme peut converger trés rapidement si m > 3, tellement qu’en pratique a 
partir d’une certaine précision il peut étre ignoré. Par exemple avec ¢(3), on a donc avec 
m = 256, on obtient une précision de 353 chiffres décimaux. 
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Les coefficients (ici C dans l’équation) est trés facile a identifier, le B est plus 
problématique mais calculable quand méme. C’est une expression polynomiale simple. La 
précision de |’approximation croit de facon géométrique, suffisant pour en extirper une 
formule pratique. 


On peut en exploiter les caractéristiques puisque : 


co 


1 713(k* + 20k2 + 16) + 360k? ¥ 1 


n3(em4/k — 1) 
n=1 


$B = — 755 k(k2 +4) 


Approxime ¢(3), ak décimales. 
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On remarque que dans le cas de k = 2, la formule est presque la méme que celle de 
Ramanujan sauf pour le facteur 2 dans la somme infinie. Pour avoir la k’iéme décimale il 


suffit alors de poser 


[10*(f(k))] mod 10 


f (k) étant l’une des formules du tableau. 


Elles convergent assez lentement mais au moins on a une formule compacte et close. 


Si on applique le méme procédé pour les autres valeurs de ¢(2n + 1), les coefficients 
deviennent assez volumineux. Rien que pour ¢ (5), je n’ai pas réussi a obtenir l’équivalent 
qui soit avec une expression polynomiale de degré 10, c’est encore pire avec ¢(7). 


Une autre possibilité est de trouver une facon d’obtenir explicitement l’expression du 


terme (€quation 1). 
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Dans un article précédent, je donnais les formes explicites des séries avec un k négatif. 
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En s’inspirant de ce modele, a l’aide du programme Pari-GP (lindep) ou PSLQ (version 
Maple) on peut identifier l’expression exacte de S(5,3) et méme S(5,a/b). En effet, pour 
S(5,3) on obtient bien (basé sur le méme modele). 
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Et de la, certaines valeurs fractionnaires 
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Ce qui donne un espoir de pouvoir calculer explicitement la valeur de S (-5, =). Ici le 
modeéle utilisé plus haut ne fonctionne pas parce que (principalement) les valeurs aux 
arguments négatifs de la fonction S (=n, “) ne sont pas simples. Par exemple S (-1, *) est 
connu explicitement : 


1 5 1 1 3 1 


Le degré du terme algébrique est trés imprévisible et les coefficients aussi. On a bien une 
x 3 ; 

expression avec (toujours les mémes) : In(2) , In(z), z et InT (2), le dernier morceau est 

inconnu et probablement de degré élevé, donc difficilement détectable. 
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